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The axio-electric effect in silicon atoms is sought for solar axions appearing owing to 
bremsstrahlung and the Compton process. Axions are detected using a Si(Li) detector placed in a 
low-background setup. As a result, a model-independent constraint on the axion-electron coupling 
constant \gAe\ < 2.2 x 10~ 10 has been obtained, which leads to the bounds tua < 7.9 eV and tua < 
1.4 keV for the mass of the axion in the DFSZ and KSVZ models, respectively (90% C.L.). 
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I. INTRODUCTION 

The axion is a neutral pseudoscalar particle introduced 
in the theory in order to explain the absence of the 
CP violation in strong interactions Q. The effec- 
tive axion-photon (gAj), axion- lepton (gAe), and axion- 
hadron (qan) coupling constants, as well as the mass of 
the axion uia, are inversely proportional to the energy /a 
at which Upq{\) symmetry introduced in [l[ is broken. 
The initial model of the standard PQWW -axion, where 
f a = {VZGf)^ 1 ^ 2 is fixed at the electroweak scale, was 
reliably excluded by reactor and accelerator experiments 
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New models of an invisible axion with an arbitrary 
scale of Upq symmetry breaking stimulated continuation 
of the experimental search for a pseudo-scalar particle 
with a mass in a wide range from 10 -16 to 10 6 eV. The 
models of the KSVZ axion (hadron axion) [f| i] and the 
DFSZ axion (GUT axion) [7J, |8| are the two main models 
of the invisible axion. The mass of the axion tua in both 
models is expressed in terms of the properties of the ir° 
meson as 



m A = 



Sa 



(l + z + w)(l + z) 



1/2 



(1) 



where m w and f v are the mass and decay constant of 
the 7T° meson, respectively, and z = m u /md « 0.56 and 
w = m u /m s w 0.029 are the quark mass ratios. 

The most stringent bounds on the interaction of the 
axion with matter were obtained from astrophysical and 
cosmological data. Astrophysical constraints based on 
the interaction of axions with photons and electrons lead 
to the bound tua < 10~ 2 eV if the standard relation be- 
tween gAe, gA-y, and tua is used Hi]. The SN1987A data 
provide upper and lower bounds for the axion-nucleon 
coupling constant [T(|- These bounds allow the axion 
mass range ps 10 eV (the so-called hadron axion win- 
dow). At the same time, the latest cosmological data 
imply that tua < 1 eV 

Experimental bounds on the mass of the axion fol- 
low from constraints on the gAy, gAe-, and gAN coupling 
constants, which significantly depend on the theoretical 
model used. In particular, the hadron axion does not 
interact with leptons and usual quarks at the tree level. 
This leads to the strong suppression of the interaction 



of the axion with leptons, which is due only to radiative 
corrections. In some models, the axion-photon coupling 
constant gA-y is strongly different from the initial values 
for the models of the DFSZ and KVSZ axions 0- 

This work is devoted to the search for axions appear- 
ing in the Sun owing to the axion-electron interaction, 
namely, bremsstrahlung from electrons in the field of a 
nucleus, e~ + Z — > Z + A, and the Compton process, 
7 + e~ — > eT + A [T^-[2Qj ■ The reaction cross sections 
and, therefore, axion fluxes are proportional to g\ e . 

To detect axions, we used the axio-electric absorption 
of an axion by atoms, which is an analog of the photo- 
electric effect. Since the reaction cross sections are pro- 
portional to g\ e , the expected count rate of axions de- 
pends only on the axion-electron coupling constant and 
is proportional to g Ae . The dimensionless constant gAe 
is generally related to the electron mass and /a as gAe — 
C e m e / j a, where C e is the model -dependent parameter 
of about 1. In the DFSZ axion model, gAe is related to 
the electron mass m e as 



gAe 



(l/3)(cos 2 /3)m e // / 



(2) 



where j3 is an arbitrary angle. In the KSVZ axion model, 
gAe is much smaller (by a factor of about a 2 ), because it 
is determined only by radiative corrections [l9j: 



gAe 



3a 2 nm e 



E f A 24 + z + w, A 
— In In — 

N rrie 3 1 + z + w m e 



(3) 



where n is the number of generations, E/N is the model 
dependent ratio of the order of unity and A pa 1 GeV is 
the QCD cutoff scale. 



II. SPECTRA OF SOLAR AXIONS AND THE 
CROSS SECTION FOR THE AXIO-ELECTRIC 
EFFECT 

The spectrum of Compton axions (c?<I>^ omp / '(IEa) ap- 
pearing in the 7 + e~— > e~ + A reaction was calculated 
by integratin g th e Compton conversion cross section ob- 
tained in (2TI |22| over the Planck energy distribution of 
thermal photons taking into account the radial distri- 
bution of the temperature and electron density in the 
standard solar model BSB05 23] for the high-metallicity 
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FIG. 1. Spectra of solar axions produced by the 
bremsstrahlung and the Compton process as calculated for 
various values of the mass of the axion and constant gAe = 
lO" 10 . 



case [24[. The calculation results for various masses of 
the axions are shown in Fig[TJ 

The spectrum of bremsstrahlung axions 
(d^ A h ~ em / '(LEa) was calculated by the same method. 
The cross section for the e - + Z — > Z + A reaction was 
obtained in (25|. The Maxwell-Boltzmann distribution 
was used to describe the energy spectrum of electrons. 
Bremsstrahlung was considered in the field of H, 4 He, 
3 He, 12 C, 14 N, and 16 nuclei. The spectrum of 
bremsstrahlung axions is softer (FigfT|). 

The calculation procedures were described in detail in 
[26j , where simple parameterizations describing the spec- 
tra with an accuracy of 1% in the range of (1-10) keV 
were presented. The spectra of axions shown in Fig. 1 
were calculated for gAe = 1CT 10 . The average energies of 
axions are 1.6 keV and 5.1 keV for bremsstrahlung and 
Compton axions, respectively, for the case rriA ~ 0. The 
maximum intensity is reached at energies of 0.6 keV and 
3.8 keV, respectively. The axion flux almost vanishes at 
energies above 20 keV. 

An axion interacting with an electron should undergo 
axio-electric absorption, which is an analog of the photo- 
electric effect. Silicon atoms entering into the composi- 
tion of a Si(Li) detector were used in our experiment as 
targets for the axio-electric effect. The cross section for 
the axio-electric effect for nonrelativistic axions is pro- 



portional to the cross section for the photoelectric effect 
a pe for photons with the energy equal to the mass of the 
axion 1211: 



a ae (E A ) = a pe {E 1 = m A ) 



9Ae 



3m 2 A 



f3 I6iram2 



(4) 



where /3 = v/c = pa/ ' Ea is the velocity of the axion. For 
rclativistic axions in the case Ea < rn e and vtla — > 0, the 
cross section differs from Eq. (4) by a factor of about 2/3 
and by a change of itia to Ea [2l|, [27J . The cross section 
for the photoelectric effect cr pe (E) for various energies of 
photons was calculated using the XCOM database [28j . 
Since the efficiency of the detection of an appearing elec- 
tron and accompanying X rays by a Si (Li) detector is 
= 100% the expected spectrum of the detected energy is 
given by the expression 



dN (TP \ 

— — = a ae {EA,m.A) 
aE 



dE A 



dE A 



(5) 



To compare with the experimentally measured spectrum, 
Eq. (5) should be additionally averaged using the re- 
sponse function of the Si(Li) detector, which is well de- 
scribed by a Gaussian for photons and electrons in the en- 
ergy range of (0-20) keV. The cross section for the axio- 
electric effect calculated for itia — is shown in Fig|2] 
along with the expected spectrum from axions S(E, mA) 
with and without the inclusion of the energy resolution 
of the Si (Li) detector (a = 115 eV). The spectrum has a 
characteristic feature at an energy of 1.84 keV, which cor- 
responds to the binding energy of electrons in the K shell. 
The peak shape of the feature significantly increases the 
sensitivity of the experiment compared to the case of a 
smooth curve. 

The number of detected electrons and X-ray photons 
following the absorption of an axion depends on the num- 
ber of silicon atoms in the detector, the time of measure- 
ments, and the efficiency of the detector. The probability 
of observing a signal is determined by the background of 
the experimental setup. 



III. 



EXPERIMENTAL SETUP 



In the experiment, we used a Si (Li) detector with a 
sensitive-region diameter of 17 mm and a thickness of 2.5 
mm. The detector was placed in a vacuum cryostat with 
the input beryllium window 20 /im thick. The window 
was used for energy calibration and determination of the 
detection efficiency of gamma-ray photons in order to find 
the sensitive volume of the detector. 

The detector was surrounded by 12.5 cm of copper and 
2.5 cm of lead, which reduced the background of the de- 
tector at an energy of 14 keV by a factor of 110 as com- 
pared to the unshielded detector. In order to suppress 
the background from cosmic rays and fast neutrons, we 



3 



f10" 1 




E A , keV 



FIG. 2. Cross sections for the (1) photoelectric effect for 
a silicon atom and (2) axioelectric effect for gAe ~ 1 and 
rriA = 0, (3) the expected spectrum of events detected per 
day in 1 g of Si in 1 keV interval, and (4) the spectrum taking 
into account the resolution of the detector (right scale). 

used five scintillation detectors, which closed the detec- 
tor almost completely except for the bottom side, where 
a Dewar vessel with liquid nitrogen was placed. 

The spectrometric channel was organized as follows. 
A field-effect transistor was mounted on a Teflon plate 
near the Si(Li) detector and was cooled to a tempera- 
ture close to liquid nitrogen. A preamplifier except for 
the first cascade was placed outside the passive shielding 
in order to reduce the background. A signal after the 
preamplifier was guided to two amplifiers with different 
amplification coefficients, which made it possible to mea- 
sure the spectrum of signals in the low energy (0.5-60) 
keV and harder (10-500) keV regions. The measurement 
of two regions allowed the reliable control of the back- 
ground in the natural radioactivity region. Each spec- 
troscopic channel was equipped with its own analog-to- 
digital converter. Thus, four 4096-channel spectra were 
simultaneously accumulated in computer memory (two 
energy ranges in coincidence and anticoincidence with 
an active shielding signal). 

The sensitive volume of the detector was determined 
using 7 lines of standard 57 Co and 241 Am sources. The 
determined number of silicon atoms in a sensitive volume 
of 0.57 ± 0.03 cm 3 was N si = 2.85 x 10 22 . The error is 
mainly due to the inaccuracy of the position of the 57 Co 




FIG. 3. Spectra of the Si(Li) detector measured in (1) antico- 
incidence and (2) coincidence with the active shielding signal. 



source with respect to the detector. 



IV. RESULTS 

Measurements continued for 76.5 days of pure time in 
the form of two-hour runs in order to control the stabil- 
ity of the operation of the spectrometric channels of the 
Si(Li) detector and active shielding scintillation detec- 
tors. Figure [3] shows the total energy spectrum of signals 
in the range of (0.8-45) keV. The spectrum of events 
detected in anticoincidence with active shielding signals 
clearly exhibits the 13.9 and 17.8 keV Np X-ray lines 
attributed to decays of 241 Am. The spectrum of events 
detected in coincidence with active shielding signals has 
no statistically significant peaks. 

The spectrum measured in the range of (0.8-9) keV 
is shown in FigU along with the expected spectrum 
S(E, tua) in the case of the detection of theaxioelectric 
effect for tua ~ and gA e — 4 x 10 -10 (line 4 in Figf2]). 
It can be seen that the shown characteristic feature is 
not observed in the measured spectrum. The maximum 
likelihood method is used to determine the intensity of 
the response function. The experimental spectrum was 
represented in the form of the sum of an exponential de- 
scribing the smooth background and the response func- 
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FIG. 4. Spectrum of the signals of the Si(Li) detector with the 
optimal fit and (solid line) the expected spectrum in the case 
of the detection of axions with tua ~ and gAe — 4 x 10~ 10 . 



tion for axions S(E, m A )'- 

N(E) = a + b exp(cE) + g Ae S{E, m A )N St T (6) 

Here, N$i is the number of silicon atoms in the sensitive 
volume of the detector and T = 6.61xl0 6 is the live time 
of measurement. 

The energy scale and energy resolution a were de- 
termined from peaks manifested in the measurements 
and were fixed. The constant g\ e and the parameters 
a, b, and c that describe the continuous background 
were free parameters. The fit in the range of (1.2- 
8.0) keV for = is shown in Fig|4j The mini- 
mum value x 2 = 111/108 corresponds to the value g\ e 
= (-0.06 ± 15) x 10 40 . To determine the upper bound 
for g Ae , we used the standard method: the x 2 value was 
determined by fixed g\ e values, whereas the other pa- 
rameters were free. The resulting probability function 
P(x 2 (9Ae)) was normalized to unity for values g\ e > 0. 
The upper bound for \gAe\ is 

\g Ae \ < 2.2 x 10- 10 (7) 

at 90% C.L. Limit (7) is a model independent bound for 
the coupling constant of the axion or any other pseu- 
doscalar relativistic particle with the electron. 

If the mass of the axion is several keV, the expected 
spectrum in the detector changes significantly and de- 
pends on the particular value. To determine values, 



FIG. 5. Spectrum of the signals of the Si(Li) detector in the 
range of (1-16) keV with the optimal fit for tua ~ 5 keV. 
The energies of the Gaussian peaks are given near them. The 
expected spectrum is shown for the case of the detection of 
axions with ttla = 5 keV and gAe — 5 x 10~ 10 . 



we used the procedure similar to that described above. 
The spectra of solar axions were calculated for different 
values of m A from 1 to 10 keV with a step of 1 keV. The 
cross section of the axioelectric absorption was calculated 
by the formula 



Cabs 



(Ea) 



J pe 



9Ae 



3E A 



f3 167ram 2 



1- ~ 



(8) 



At P — > 1 and j3 — > 0, this formula coincides with the 
cross sections for relativistic and nonrelativistic axions 
obtained in [27] (Eq. (4)) and provides an extrapola- 
tion approximately linear in /3, which ensures a sufficient 
accuracy for the case under consideration. 

An example of the response function of the detector 
for detecting axions is shown in FigJS] The spectrum 
fitting range was expanded to 16 keV. To describe the 
experimental spectrum in a wide range, function (6) was 
supplemented by a linear term for describing the con- 
tinuous background and six Gaussians for describing the 
peaks of the characteristic Np X rays manifested in the 
measurements [29|, |3(| • The fit for tua = 5 ke V is shown 
in FigH 

The maximum deviation of g\ from zero for all m A 
values is 2.5a. The upper bounds obtained for \g A e\ a t 



5 





\ 4 \ 






1 5 -| 




yf 


- 

3~ 


10" 7 -! 


\ 9 / 






< 

CD 

-Q 

10- 


1 - present N/ 

work / y* 




- 




/ 7 a y 






10 11 - 


/% z / s^J 

/ KSVZ/ Y 

/ / 8 






-lO" 13 " 


I 1 '"""I 1 1 1 


i i 





10" 2 10° 10 2 10 4 10 6 

m A , eV 



FIG. 6. Bounds for the axion-electron coupling constant from 
(1) this work, (2) reactor experiments and solar axions with 
energies of 0.478 and 5.5 MeV [H[-[3i|, (3) beam dump ex- 
periments Hlli (4) d ecay of orthopositronium [3|§], (5) 
CoGeNT [3^ (6) CDMS [1|, (7) bound for the axion lumi- 
nosity of the Sun [23], (8) red giants [13|, and (9) experiment 
with 169 Tm [2(J. The regions of excluded values lie above the 
corresponding lines. The inclined lines show the gAe values 
in the DFSZ and KSVZ (E/N = 8/3) models. 



various niA values are shown in Fig [6] (line 1) in compar- 
ison with the other experimental results. 



For large masses of the axion tua > 10 keV, the most 
stringent bounds were obtained from reactor [3lL I32I ] and 
accelerator [36l.l37| e xperiments, as well as for high energy 
solar axions |33l - l35f . For small rriA values (tua < 0.1) 
keV, the most stringent bounds were obtained from astro- 
physical data, namely, from the cooling rate of red giants 
|41j and under the assumption that the axion luminosity 
of the Sun does not exceed 10% of the observed neutrino 
luminosity 22] . The search for axions with masses in the 
intermediate region wia ~ 1 keV was specially performed 
in view of the possibility of attributing the DAMA ex- 
periment result to the axio-clcctric effect with silicon and 
ger manium detectors designed for seeking dark matter 
[3^ . HqJ . Solar axions having a continuous energy spec- 
trum should lead to the resonant excitation of low-lying 
nuclear levels [26| . 

The upper limit obtained for \gAe\ lead to bounds 
for the mass of the axion that depend on a particu- 
lar model. In the model of the DFSZ axion, the rela- 
tion between gAe and tua for cos 2 /? = 1 has the form 
gAe = 2.8 x \Q~ ll rriA, where itia is expressed in eV units. 
The limit (7) leads to the bound m% FSZ < 7.9 eV (at 
90% C.L.). The upper bound for the excluded values 
units is 15 keV (Fig©. In the model of the KSVZ axion, 
gAe & n d niA are related through Eq. (3). For the ratio 
E/N = 8/3, which is characteristic of most variants of 
the DFSZ axion, and n — 3, the region of excluded values 
mf svz is (1.4-12) keV (Fig©. 

V. CONCLUSION 

The axio-electric absorption of solar axions, which pro- 
duced by the Compton process and the bremsstrahlung, 
by silicon atoms has been sought. To detect the 
axio-electric effect, a Si(Li) detector placed in a low- 
background setup was used. As a result, new limits have 
been obtained for the axion-electron coupling constant 
and mass of the axion. For axions with a mass smaller 
than 1 keV, the resulting limit is \gAe\ < 2.2 x 10~ 10 at 
90% C.L.. 
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